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ON THE INFLUENCE OF KEYWAYS ON THE STRESS DISTRIBU- 
TION IN CYLINDRICAL SHAFTS* 

BY 

T. H. GRONWALL 

1. When a keyway is cut in a cylindrical shaft, subject to torsion only, 
and of circular cross section, the maximum stress in the keyed part will con- 
siderably exceed that in the full part. To form an estimate of this increase, 
Filonf has investigated the torsion in a shaft the cross section of which is 
composed of confocal ellipses and hyperbolas. His results are expressed in 
infinite series of trigonometric and hyperbolic functions, and their numerical 
computation is necessarily somewhat laborious. Moreover, his numerical 
results (1. c, § 22-23) extend only up to a ratio of minor and major axes in 
the ellipse equal to tanh (7r/2) = 0.917, and this section differs too much 
from a circle to allow a sufficiently close numerical estimate of the stress 
ratio in the latter case. 

In the present paper, a cross section will be considered which frequently 
occurs in practice: the shaft section is taken as a circle, and the keyway is 
formed by a circular arc intersecting the former orthogonally. 

We shall first determine the exact formulas for the stress distribution, and 
then derive from these an approximate formula, adapted to practical use, in 
which the ratio of the maximum stresses in the keyed and full parts of the 
shaft is represented by the expression 
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where a and b are the radii of shaft and keyway respectively. This expression 
is exact for 6 = (and thus shows that a flaw or crack in the surface of a 
circular shaft has the effect of doubling the maximum stress) but is somewhat 
too great for < b/a ^ j, this range of the ratio b/a being sufficient for all 
cases occurring in common practice. 

* Presented to the Society, Feb. 26, 1916. 

t L. N. G. Filon, On the resistance to torsion of certain forms of shafting, vdth special 
reference to the effect of keyways, Philosophical Transactions of the 
Royal Society, London, ser. A, vol. 193 (1900), pp. 309-352. 

234 



STRESS DISTBIBUTION IN CYLINDRICAL SHAFTS 



235 



If a linear formula for the ratio of maximum stresses in the keyed and full 
parts of the shaft is desired, the above expression may be replaced by 

^ 5a' 

which is greater than the preceding one for < a/6 ^ J. 
2. With the notations of the figure, it is seen at once that 



(1) 



tan a = - , 
a 




Fig. 1. 



and that the distance between centers of shaft and keyway is a sec a . 

The first step toward the determination of the torsional stresses consists 
in the construction of a harmonic function ^ ( a; , y)* taking the value 
^{x^ + y^) on the boundary of our cross section. To this purpose, write 
z = x + iy,w = u + iv, and map the area ABCD conformally on the first 
quadrant in the w-plane in such a manner that C corresponds to w = , D to 
w = <x> , the circular arc CBD to the w-axis, and the arc CAD to the i)-axis. 
It is seen immediately that 



(2) 



z — ae 



w = 



e""z 



z = a 



w — e' 
'w - 1' 



dw 
dz 



1 {e'^'w - ly 
a e^' - 1 ' 



and the segment BA of the real axis corresponds to the first quadrant of the 
circle -w^ + «^ = 1 . From (2) we obtain 



(3) l(a;2 + 2/2) =1|2|2 = 



a^ 1 — 2 cos a-u — 2 sin a-v -\- v? -{-v^ 
2 1—2 cos a • -M + 2 sin a • » + M^ + t)^ ' 



* The notations are those of Love's Theory of elasticity, 2d edition, chapter XIV, except 
that his J , >7 , X, , Yz are replaced hy u ,v , X , Y . 
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so that the boundary conditions for ^ become 

^a^ for t) = 
(4) ^ = i\z\^ = ' 



1 z/^i isma-v \^ 

ia^ I 1 - -, . o • 1 — 2 I for M = . 

\ 1 + 2 sin «•■!) + t)''/ 



Now form Green's function G{u' , v'; m, ») in the first quadrant of the w- 
plane (u' ,v' are the coordinates of the variable point and u , « those of the pole). 
By reflecting the pole successively into the positive v-, the negative u-, and 
the negative t)-axis, we obtain at once 

, [ j u' -uy + w + «)^] [ {u' + uf + (/ - ^ )^] 

The formula 

dG 



^ 27r J ^'dn 



ds 



now gives 

22r _ r°° 4sinQ;-/ T u u "1 

(5) -^^ = 7r-j^ 1 + 2 sin «•»' + «'^ U^ + («' - s)^ ~ m^ + («' + »)^J ' 

whence, introducing the function $ conjugate to ^ , 
$ + i^ = 



(6) + 



2 

2a^ tan a tv — 2a — it cos a-w— (tt — 2a) cos 2a-w'^ -\- ir cos a • w^ 
TT 1—2 cos 2a- vP" + w* 

8a^ sin^ a w^ log w 



■K 1—2 cos 2a- V? -\- w* 



for < a < ■7r/2 , where log w is the principal value of the logarithm. 
For « = , w = M and we obtain from (6) 

2a^tana:r cosa-M— 1 cosa-M + 1 "I 

$ = a z 7. ; -, + ( TT — a ) :; — —z , i 

TT L 1 — 2 cos a-M + w' 1 + 2 cos a-M + w'J 

(7) 8a^ sin^ a m^ log u 

'v 1 - 2cos2a-w2 + M*' 

T «' 
^ 2 ' 

while f or ■M = , w = zV , log w = log » + 7ri/2 and 

_ 2a^ tan a (y — 2q:) (1 + cos2q;-2)^) 8a^ sin^ a v^ log » 

^ ~ Z 1 i o r>.. ..2 i TJ r 



TT 1 +2cos2o;-z)2 + »* ' tt 1 + 2 cos2a-»2 + 1)«' 

(8) 

a? 2a^sina:-» 

2 ~ 1 + 2sina-» + 2)2 ' 
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the identity of these boundary values of ^ with those given in (4) serves as a 
check on the calculation. 

3. We now proceed to determine the maximum stress in the cross section. 
The stress components X and Y are given by 

Y = ^r(py+ x) = - ^r(^~ - xy 

and the maximum stress is the maximum of VZ^ + Y^ , which by an elementary 
property of harmonic functions must occur on the boundary.* For m = 0, 
we have 

and since "^ — ^ {x^ -\- y^) =0 for m = 0, the partial derivative of this 
expression in respect to v vanishes, so that 



Similarly 



whence 






From (8), (3), and the last of (2), we now obtain 



(9) 



jura 

_ 2(7r - 2«)t)(cos2a + 2t)^ + cos2a-t)^) 
(1 + 2 sin a • t) + «^ ) ( 1 — 2 sm. a-v -\- v^Y 

2sin2a;[t)(l + 2 cos 2a; • t)^ + t;^ ) +2^(1 - v'^)\ogv] 
(1 + 2 sin a • ■» + t)^ ) ( 1 — 2sma-xi -\- v^Y 

TT ( 1 + cos 2a) V 



+ 



1 + 2 sin a • ■» + t)^ ■ 



' See Love, Theory of elasticity, 2d ed., § 219. 
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The stress distribution being evidently symmetrical about the a;-axis, the 
stresses on the arc AC are, in particular, equal to those on the arc AD, and 
this circumstance finds its expression in the fact that (8) is invariant when v 
is replaced by its reciprocal. 

From /(») =/(!/■») it follows that /'(f) = — f {l/v)/iP, whence 
/' ( 1 ) = , so that fiij) becomes a maximum or minimum at » = 1 . To 
find any other maxima or minima, it is sufficient to consider the interval 
^ t) ^ 1 ; writing 

so that t ^ 1, (9) takes the form 

Ef - 2Ft + G + 2sm2a-t 



2/(») = 



(® --yogv 



{t + sin a) {t — sin ay ' 

where 

£ = TT + ( Sir — 4a ) cos 2a — 2 sin 2a; , 

F = IT sin a ( 1 + cos 2a ) , 

G = sin^ a ( Stt — 4a + IT cos 2a + 2 sin 2a ) . 
The equation/' (v) = now becomes 
[2(f -sin^a)[{E + 2sin2a)t- F] 

- {St + sin a)[Ef - 2Ft + G]} (^1 -l,^ 

+ 2 sin 2a-log» { {f - sin^ «)[(«" ^) (l " ^2) 

or reducing by (10) and dividing by 1 — 1/t)^, thus removing the known 

root V = 1 , 

X(«) =2{fi - sin^ a)[{E + 2 sin 2a) t- F] 



(11) 



- idt + sin a)[Ef - 2Ft + G ] 

+ 8 sin 2a [ {fi - sin^ a)i2f -1) 

loe v 

- (St + sin a) {f - t)]-^ = . 



V 

We may prove that this equation has no root between and 1 as follows. 
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Since any root of x ( ^ ) =0 in the interval < ?> < 1 varies continuously 
with a , three cases are to be considered when we let a decrease steadily from 
the given value to zero : (a) there exists some value ao S of a such that the 
root in question approaches zero as a -» ao; (&) the root in question approaches 
the value v = \ for some value of a , or (c) the root will remain within the inter- 
val < « < 1 as a decreases toward zero. It is readily shown that all three 
cases are impossible. In case (a) we find from (11) 

1 

limsup^^s^ = [2(E + 2sin2a) - 3E ],=„„ - 4 lim inf sin 2a • ^- log-, 

«._»0 " v-^0 1 V 

a'^ao a— >ao — — V 

V 

and since 

1 
-+^1,1 



log - > log - > 1 for V < e ^ , 



we have 



1 V ^ V 

V 



1 

■V 1 

lim inf sin 2a- z log- S sin 2ao, 



SO that finally 

X (^) ■"" 

lim sup — ij— m - Ea=ao < for ^ oro < o • 

In case (b), we obtain at once 

x(l) = — 2cos2a [(2 + 4 sin a) (ir — a) — 4 sin 2a 

+ sin^ a ( 9 — 2 cos a — ir sin a — 8a ) ] , 

and since 9 — 2 cos a — ir sin a — 8a decreases as a increases from to ir/2 , 
we have 

X ( 1 ) < — 2 cos^ a [ ( 2 + 4 sin a)(ir — a) — 4 sin a— (Sir — 9) sin^ a ] . 

In the interval ^ a = ir/G , we find 

(2 + 4sina)(ir - a) - 4 sin a - (Sir - 9)sin2a > 2(ir - ir/6) 

-4-i- (5ir-9)-l>0, 

and a similar argument applies to each of the intervals ir/6 = a = ir/3, 
ir/3 ^ a ^ arcsin 4/S, arcsin 4/S ^ a ^ ir/2, so that finally 

X(l) <0, O^a < ir/2. 
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Hence /" ( 1 ) < , so that » = 1 gives a maximum of ■yX^ + Y^ . In the 
third case, it is seen that for < t; < 1 

limx(«) = - 4:Tr^ < 0, 

a-»0 

SO that none of the three cases yields a root of (11). For » = 0, (9) shows 
that VX^ + y^ = 0, and since x ( ^ ) =0 has no root in the interval < « < 1 , 
the maximum of VX^ + Y^ occurs at t) = 1 for ^ a < 7r/2 , and its value is 
found from (9) to be 

a2^ M J V2 I y2 = 1^ 2 {it - a) - 27r sin a + tt sin^ « - 2 sin « cos « 

TT COS a ( 1 — sin a ) 

Referring to the figure, the maximum stress on the arc CAD is therefore seen 
to occur at A , while the stress is zero at C and D . On the remaining part of 
the boundary, v = , and here we find, in the same way as for m = , 

Treating this expression in exactly the same way as the preceding one, we find 
a unique maximum value of the stress at w = 1 (corresponding to the point 
B ) , which however is smaller than (12), the latter therefore giving the absolute 
maximum. 

4. To determine the constant t, we have the following relation, where M 
is the moment of the external forces about the axis of the shaft, 

M = 2nT j j{<b - i(x' + y')]dxdy. 

Denoting by 7(m, v) the integral in (5), this becomes 

M = fira^ I I dxdy — fiT I I {x^ + y^) dxdy I I I{u,v)dxdy , 

or changing the integration variables in the last integral and remembering 
that 

d{x,y) dz 

d {u, v) dw 

which may be found from (2), we obtain 

M = ixra^ I I dxdy — fXT i | (x^ + y^)dxdy 

(13) 

ifXTO^ sin^ a f"" f"" I{u,v)dudv 

TT Jo Jo (1 — 2 COS a-M + 2 sin «•■!) + w^ + j)^)^ 
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The first integral in (13) equals the area of the cross section, and is most 
readily found by drawing CD and evaluating the area as the difference of 
the two circular segments CDBC and CDAC, of radii a and a tan a and 
center angles 27r — 2a and w — 2a respectively. Thus we find 



(14) 



//^^,= (,-„ + si„<.cos„,„>-(|-„-oo.<. sin <.)«'*.„', 

= a^ ( TT — a + tan a — ^ tan^ a + a tan^ a I . 

The second integral in (13) is the polar moment of inertia; introducing polar 
coordinates r,d, the equation of the circular arc of radius b is 



r = a sec a ( cos 6 — Vcos^ 6 — cos^ a ) , 
and we have 

f f{x^+y^)dxdy = 2 f dd f f" dr + 2 f dd f r^ dr 

= 5a^(7r-a) + 5a'*sec^a I ( cos ^ - Vcos^ 6 - cos^ afdO 

Jo 

= |a*(7r- a) + |a*sec^a f ( 8 cos^ 9 - 8 cos^ cos a^ ^ 

Jo 

+ COS* a)de + 2a^ sec^ a j cos ^ (2 sin^ 6-1 

— sin^ a ) Vsin^ a — sin^ 9 dd . 

Making sin ^ = sin a sin <j> in the last integral, the integrations may be 
performed, and expressing all trigonometric functions in tan a , we obtain 



+ a tan^ a — ^ tan^ a + 3a tan* 



«]■ 



In the last integral in (13), I{u,v) may be obtained explicitly from (6), but 
the further integration in respect to u and v is not possible in finite terms. 

5. We shall now derive approximate formulas, adapted to numerical calcu- 
lation, from the purely theoretical results of the preceding paragraphs, and 
begin with the expression for the maximum stress. Equation (12) may be 
written 

Max. ^^X^TY' = 2y.Ta - ^^ tan a + ^"^ 



X 



TV TT 

[27r — (27r — 2)sina](l — COS a) — 2a + 2 sina + (tt — 4)sin2a 

cos a ( 1 — sin a ) ' 
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replacing a by sin a , the last expression is increased (except for a = ) and 
we find 

4jura 



Max. VX2 + F2 ^ 2ina - -^^tan a 



IT 




1 ''^"'t;,,.' r(27r-2)cosa 


2 cos a 


1- ^ tan ..|^ i+eosa ' (1- 


- sin a ) ( 1 + cos a ) 


( TT — 4 ) cos a " 
1 — sin a J ■ 





Now each of the three terms in the bracket decreases as a increases (as is 
seen upon forming their derivatives), and the value of the bracket for a = 
being 27r — 4 , we obtain 

(16) Max. VZM^^ ^ MTffl [2 - ^tan a + ^2 - ^ jtan^a I . 

Next, we proceed to find an upper bound for the third integral in (13). From 
(5) it follows that, the expression in brackets under the integral sign being 
evidently positive, we have for < a < 7r/2 

r, , . r 4^' r ^ ^ 1^, 

I(u, V) < sma I -i\ o I / / ^ — ^rn — rr^, — \2 "^ ' 

^ Jo I + v' lu^ + (v' - vy u^+iv' + vy} 

and performing the integration, 

4irsina;''» 
( 1 + M )^ + r 
Hence 

r" I{u,v)dv 

J, (1 — 2cosa-M + 2sina-» + M^ + 2^)^ 

Air sin a-vdv 



^Jo for 



TT sin a r 1 1_ 

~ 1 + cos a\_u 1—2 cos a-u ■{■ v? 

1 1 1 — 2 COS«-M + M^ "| 

+ 2~(l + cosa) 'u^^^ (1 +m)2 J 

and furthermore 

r" r* I{u,v)dudv 

Jo Jo (1 — 2 cos a-M + 2 sina-'i) + M^ + «^) 

TT sin g r°° r 1 1 1 

1 + COS a Jo \_u 1 — 2 cos a-u + u^ 2(1 + cos a ) 

1 1—2 cos a-u + u^~\ 
._log ___ ^du. 
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Integrating by parts so as to remove the logarithm, the last integral is 
readily evaluated, and we finally obtain 



/loo /!« 

Jo Jo (i 



I {u, v)dudv 



2 cos a-u + 2 sin a-v + u^ + v^y 
(17) ----- 

7r(x — a — sina;)(l — cos a ) 



for = a < 7r/2, the equality sign holding in (16) and (17) for a = only. 
From (13), (14), (15), and (17) it is now seen that 

M ^ —^ [w — a + tan a + a tan^ a — 3 tan^ a + Sir tan^ a — 3a tan^ a 

— 8 ( TT — a — sin a ) ( 1 — cos a ) ] . 
We now use the inequalities 

tan a — I tan^ a < a < tan a — ^ tan^ « + i tan^ a < tan a , 

tan a , , 

sin a = — p — > tan a — 5 tan'* a , 

Vl+tan^a 

1 — cos a = 1 , — = < i tan^ a , 

Vl+tan^a 



valid for < tan a < 1, and obtain 



fxra* 



M S —^ [ T — ( tan a — I tan^ « + i tan^ a ) + tan a 

+ ( tan a — J tan^ a ) tan^ a — 3 tan^ a + 37r tan^ a — 3 tan a • tan* a 
— 8 ( TT — ( tan a — i tan^ a) — { tan a — | tan' a ) ) • 5 tan^ a 



fiTTa* 



\ 1 — tan^ a 4 — ^ tan a ( 1 + Tq tan a — -j^ tan^ " ) | 



Introducing the assumption ^ tan a = }, whence S a ^ 14° 02' 10".48, 
the innermost parenthesis is positive, so that finally 

(18) MS^^(1 -4tan2a). 

From (16) and (18) it is seen that, replacing tan a by b/a , for ^ b/a ^ J , 

2M ^ ira^y w J a^ 



Max. VZ2 + P ^ 



T«' , . b' 

1-4^ 
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For the full part of the shaft (cross section a circle of radius a ) we have the 
well known formulas 

TV or 2 

and consequently the ratio of maximum stresses in the keyed and full parts of 
the shaft is less than 

(19) — -^^^-V^' 

for < b/a ^ J , but equal to this expression ( = 2 ) for 6 = . For practical 
use, the values of (19) may be tabulated, or we may replace (19) by a linear 
expression 2 + X ( b/a ) , where X is to be given the least value for which this 
expression is greater than or equal to (19) throughout the interval ^ b/a ^ j . 
This leads to the condition 

(<i-i)^-i 

^= b^ 

1 -4-2 
and the right hand member decreasing as b/a increases from to J , we find 



(e-i).l-i 

\ IT ) 4: IT „^,„^ 

-^ r-^^-j = - 0.2122, 

J- A 



4 

SO that the ratio of maximum stresses never exceeds 

(20) 2 - 0.2122 - 

in the interval ^ b/a ^ | . 



